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Abstract 



We study the evolution of cosmological perturbations on large scales, up to second order, 
for a perfect fluid with generic equation of state. Taking advantage of super-horizon con- 
servation laws, it is possible to follow the evolution of the non-Gaussianity of perturbations 
through the different stages after inflation. We find that a large non-linearity is generated 
by the gravitational dynamics from the original inflationary quantum fluctuations. This 
leads to a signiflcant enhancement of the tiny intrinsic non-Gaussianity produced during 
inflation in single-fleld slow-roll models. 
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1 Introduction 



Inflation is the simplest and most successful mechanism proposed to date for the causal 
generation of primordial cosmological perturbations on cosmologically relevant scales 
[1]. The gravitational amplification of the primordial perturbations is supposed to seed 
structure formation in the Universe and produce Cosmic Microwave Background (CMB) 
anisotropics in agreement with observational data [2]. Due to the smallness of the primor- 
dial cosmological perturbations, their generation and evolution have usually been studied 
within linear theory [3] . After the seminal work by Tomita [4] , only recently second-order 
perturbation theory [5-7] has been employed to evaluate specific physical observables gen- 
erated during infiation [8,9]. The importance of an accurate determination of higher-order 
statistics as the bispectrum comes from the fact that they allow to search for the signa- 
ture of non-Gaussianity in the primordial perturbations which is usually parametrized by 
a dimensionless non-linear parameter /nl- Indeed, a number of present and future CMB 
experiments, such as WMAP [10] and Planck, have enough resolution to either constrain 
or detect non-Gaussianity of CMB anisotropy data with high precision [11]. 

The main result of the second-order analysis performed in [8, 9] is that single-field 
slow-roll models of infiation give rise to a level of intrinsic non-Gaussianity which - at the 
end of the infiationary stage - is tiny, being first-order in the slow-roll parameters. 

The goal of this paper is to study the post-infiationary evolution on super-horizon 
scales of the primordial non-linearity in the cosmological perturbations. We perform a 
fully relativistic analysis of the dynamics of second-order perturbations for a perfect fiuid 
with generic equation of state taking advantage of the super-horizon conservation of the 
second-order gauge-invariant curvature perturbation recently discussed in Refs. [12, 13] 
(see also [8, 14, 15]). Our main result is that the post-inflationary evolution gives rise 
to an enhancement of the level of non-Gaussianity on super-horizon scales. Once again, 
infiation provides the key generating mechanism to produce super-horizon seeds, which 
are later amplified by gravity. 

The plan of the paper is as follows. In Section 2 we provide the second-order expansion 
of the metric and of the energy-momentum tensor, assuming that the source term is 
represented by a perfect fluid with constant equation of state. In Section 3 we solve the 
perturbed Einstein equations up to flrst order around a Friedmann-Robertson-Walkcr 
background. The body of the paper is contained in Section 4, where we derive the 
super-horizon evolution equations of the second-order gravitational potential and density 
perturbations. Section 5 contains a brief discussion of our findings. 

2 Perturbations of a flat Robertson- Walker Uni- 
verse UP to second order 

In order to study the perturbed Einstein equations, we first write down the perturbations 
on a spatially fiat Robertson- Walker background following the formahsm of Refs. [5,6]. 
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We shall first consider the fluctuations of the metric, and then the fluctuations of the 
energy-momentum tensor. Hereafter greek indices run from to 3, while latin indices 
label the spatial coordinates from 1 to 3. If not otherwise specified we will work with 
conformal time r, and a prime will stand for a derivative with respect to r. 



2.1 The metric tensor 

The components of a perturbed spatially fiat Robert son- Walker metric can be written as 
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where the scale factor a is a function of the conformal time r. The standard splitting of 
the perturbations into scalar, transverse {i.e divergence-free) vector parts, and transverse 
trace-free tensor parts with respect to the 3-dimensional space with metric 5ij can be 
performed in the following way: 
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where (r) = (1), (2) stand for the order of the perturbations, Ui and Xi ^'^'^ transverse 
vectors {d^u^^^ = d^xt^ = 0), xtj' is a symmetric transverse and trace-free tensor {d''Xij' = 
0, X^t^ = 0) and = didj — (1/3) 5ijd^dk is a trace-free operator. Here and in the 
following latin indices are raised and lowered using 5^^ and respectively. 
For our purposes the metric in Eq. ()2.H) can be simplified. In fact, first-order vector 
perturbations are zero; moreover, the tensor part gives a negligible contribution to second- 
order perturbations. Thus, in the following we can neglect ujf'\ X^(i)i a^id However 
the same is not true for the second order perturbations. In the second-order theory the 
second-order vector and tensor contributions can be generated by the first-order scalar 
perturbations even if they are initially zero [6]. Thus we have to take them into account 
and we shall use the metric 
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The controvariant metric tensor is obtained by requiring (up to second order) that g^^g^^ 
5^ ^ and it is given by 

= _Q-2(^^ _ 20(1) _ 0(2) + 4 (0(1))2 _ ^.^(D^.^d)^ ^ 

di^{^) + ^ (ai^(2) + ^i(2)) + 2 (V'(i) - 0(1)) 9^a;(i) - 9^a;(i)L>^ ^x^'^ 
+ 2V'(i) + V'^'^ + 4 (V'^^y ) - (x^'^ + ^X' 



g'' = a-^(r 



.(2) 



_ 1 (a^;^^(2) + ^•^^(2) + ^»i(2)) _ a*^(l)^u;(l) 

- 4V'(i)£>^^x(i) + L'^'^x^^^^^feX^^^] ■ (2.5) 

Using ^f^,^ and g*^^ one can calculate the connection coefficients and the Einstein tensor 
components up to second order in the metric fluctuations. They are given in the Appendix 
A of Ref. [8]. Prom now one, we will adopt the Poisson gauge [16] which is defined by 
the condition a; = x = Xi = 0- Then, one scalar degree of freedom is eliminated from 
goi and one scalar and two vector degrees of freedom from gij. This gauge generalizes the 
so-called longitudinal gauge to include vector and tensor modes and contains a solenoidal 
vector uj^'^\ 

2.2 Energy-momentum tensor of the fluid 

Since after inflation and reheating the Universe enters a radiation-dominated phase and, 
subsequently, a matter- and dark energy-dominated phases, we shall consider a generic 
fluid characterized by an energy density p and pressure P with energy-momentum tensor 

r^ = (p + p)M^M, + pn, (2.6) 

where is the four- velocity vector subject to the constraint g^^u^Uy = — 1. At second 
order of perturbation theory it can be decomposed as 

u^=\ {5^, + + . (2.7) 

For the first- and second-order perturbations, we get 
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Similarly, we obtain 
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The energy density p can be split into a homogeneous background Po(t) and a perturbation 
5p(r, x*) as follows 

p(r, x') = po(r) + 5p{T, x') = po(r) + S'^'^ p{t, x') + ^S^'^ p{t, x') , (2.10) 

where the perturbation has been expanded into a first and a second-order part, respec- 
tively. The same decomposition can be adopted for the pressure P. 

Using the expression ()2.1()j) into Eq. ()2.6|1 and calculating T'^ up to second order we 
find 

T't = Tt^°) + 5(^)rt + (J^^^Tt , (2.11) 
where T'^j^'^ corresponds to the background value, and 

Tr+5«n = -Po-5^'v, (2.12) 

5(2)n = -^5^'V-(l + «^)Pot;f^^^a), (2.13) 
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In the previous expressions we have made the assumption that the pressure P can be 
expressed in terms of the energy density as P = wp with constant w. 

3 Basic first-order Einstein equations on large- 
scales 

Our starting point are the perturbed Einstein equations ^G^^ = ^T'^ in the Poisson 
gauge. Here = SttG^. At first-order, the (0 — 0)- and {i — 0)-components of Einstein 
equations are 



-K^ (1 + w)pot;(i) 
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where we have indicated by 7i = ^ the Hubble rate in conformal time. These equations, 
together with the non-diagonal part of the {i — j)-component of Einstein equations, give 
^(1) = 0(1) and, on super-horizon scales, 

2 p„ 2 p; ' ' 

The continuity equation yields an evolution equation for the large-scale energy density 
perturbation 

+ 3n{l + w) - (1 + w)po = l^V' (^(1)' + n^P^'^) . (3.4) 

o /~L 

This equation, together with the the background continuity equation p'f^ + 3T-C (1 + w) po = 
0, leads to the conservation on large-scales of the first-order gauge- invariant curvature 
perturbation [3] 

^(1) = _ . (3.5) 

Po 

Indeed, both the density perturbation, Sp and the curvature perturbation, ip, are in 
general gauge- dependent. Specifically, they depend upon the chosen time-slicing in an 
inhomogeneous universe. The curvature perturbation on fixed time hypersurfaces is a 
gauge-dependent quantity: after an arbitrary linear coordinate transformation at first- 
order, t ^ t + 5t, it transforms as -ip^^^ — > tp^^^ + T-Cdt. For a scalar quantity, such as 
the energy density, the corresponding transformation is 5p^^^ Sp^^^ — p'^St. However 
the gauge- invariant combination (^(^) can be constructed which describes the density per- 
turbation on uniform curvature slices or, equivalently the curvature of uniform density 
hypersurfaces. On large scales C*-^'*' — 0. Using Eq. ()3.3|) and the background continuity 
equation, we can determine 

^^,„_3£+^ (3.6) 

5 + 3w 

which is useful to relate the curvature ■j/''-^-* during either the matter or the radiation 
epoch to the gauge-invariant curvature perturbation C,^^^ at the end of the inflationary 
stage. Indeed, since C*-^"* is constant, we can write 

.-^-^d", (3,7) 
where the subscript "/" means that C'^^'' is computed during the inflationary stage. 
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4 Basic second-order Einstein equations on large- 
scales AND NON-GaUSSIANITY 



In order to determine the non-Gaussianity of the cosmological perturbations after infla- 
tion, we have to derive the behaviour on large-scales of the metric and the energy den- 
sity perturbations at second order. Again, our starting point are the Einstein equations 
perturbed at second order = 5^^-'T'^ in the Poisson gauge. The second-order 

expression for the Einstein tensor can be found in any gauge in the Appendix A 

of Ref. [8] and we do not report it here. 

• The (0 — 0)-component of Einstein equations (see Eq. (A. 39) in Ref. [8]) leads to 



-87/>(i)v27/>(i) - 3 (^(1)')' = «:2a25(2)T0, . (4.1) 

• A relation between the gravitational potentials at second-order ip^"^^ and 0^^-' can be 
obtained from the traceless part of the {i-j) components of Einstain's equations (see 
Eqs. (A.42) and (A.43) in Ref. [8]). We find 



V'(2)_0(2) = _4(^(i))2_v-2(2aV«9,^«+3(l + w)7^=^t;ii)i;(i),) 

+ 3V-*^^^^ (29V(^)9jV^(^) +3{l + w) 'H^v\^f(^i)j) . (4.2) 

This constraint is the second-order equivalent of the linear constraint i})^^^ = 0'-^^ in 
the Poisson gauge. 

• In order to close the system and fully determine the variables 'j/'*-^-', 0^^'' and 5^'^^ p, 
we use the energy conservation at second-order and the divergence of the {i — 0)- 
component of Einstein equations (see Eq. (A. 40) in Ref. [8])^ 

^('V' + 37^ (1 + w;) 5(2)p _ 3 (1 + y;) Pp^(2)' _ 6(1 + w)V^«' [(5«p + 2poV^^'^] 
= -2(1 + w)po (vi'^vl^J - 2(1 +w){l- 3w)npovi'^vl) 

+ 4(1 + w)pod,tlj'-'^vl) + 2^ {^^V'^'-'^' - ^(i^'V^V'^'^) • (4.3) 

■^Notice that Eq. 14. 3() generalizes Eq. (5.33) of Ref. [13] and corrects a sign misprint in front of the 
fourth term of that equation. 
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This equation can be rewritten in a more suitable form 
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where the argument on the left-hand side can be further simplified to 

Po V Po / 
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and the final form has been obtained employing Eq. ()3.7j) [17]. From Eqs. ()4.4|1 and 
(P7S|) we find 
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where C is a constant in time, C = 0, on large-scales, and 
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4.1 Determination of the non-linearity parameter 

Since we are interested in the determination of the non-linear parameter /^^ after the 
inflationary stage, it is convenient to fix the constant C by matching the conserved quantity 
at the end of inflation (r = tA 



Pn T ^ 



Po/ 
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where we have used the fact that during inflation wj ~ — 1. 



The inflationary quantity [ipf^ + Tii^y^ ) has been computed in Refs. [8,9] 
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in terms of the slow-roll parameters e = 1 — lUi / Hj and r] = 1 + e — Hjip') where Hj is 
the Hubble parameter during inflation and (f is the inflaton field driving the exponential 
growth of the scale factor during inflation [1] . Since during inflation the slow- roll param- 
eters are tiny, we can safely disregard the intrinsically second-order terms originated from 
the inflationary epoch. 

Combining Eqs. (jHH), ()4.2p . ()4.6p and ()4.8|) we single out an equation for the gravita- 
tional potential (f)^"^^ on large scales 



= (5 + 3w)n (V^('))' + ^Hil + w) (29V(^)9,^(^) 
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where Si denotes the R.H.S. of Eq (Q. 

We want to integrate this equation from tj to a time r in the matter-dominated epoch. 
The general solution is given by the solution of the homogeneous equation plus a particular 
solution 
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where 6(r) stands for the source term in the R.H.S of Eq. ()4.10|) . Notice that the 
homogeneous solution during both the radiation and the matter-dominated epoch de- 
creases in time. Therefore we can neglect the homogeneous solution and focus on the 
contributions from the source term 6(t). At a time r in the matter-dominated epoch 
exp[— dr'TC (5 + 3w)/2] oc r~^. Thus if we are interested in the gravitational poten- 
tial (j)^"^^ during the matter dominated epoch the contributions in the particular solution 
coming from the radiation-dominated epoch can be considered negligible. Recalling that 
during the matter-dominated epoch the linear gravitational potential ip^^^ is constant in 
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time, it turns out that 
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where we have used Eq. fj3.2|) to express the first-order velocities in terms of the gravita- 
tional potential, and we have taken into account that during the matter-dominated epoch 
S[ =0. 

As the gravitational potential ijj^^^ on super-horizon scales is generated during inflation, 
it is clear that the origin of the non-linearity traces back to the inflationary quantum 
fluctuations. 

The total curvature perturbation will then have a non-Gaussian (x2)-component. For 



instance, the lapse function 

1 



can be expressed in momentum space as 
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d% d^k2 S^^^ (ki + k2 - k) /4 (ki, ks) 0(^)(ki)0(')(k2) , 
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where we have defined an effective "momentum-dependent" non-linearity parameter f^^. 
Here the linear lapse function (p^^^ = -ip^^^ is a Gaussian random field. The gravitational 
potential bispectrum reads 



(0(ki)0(k2)0(k3)) = (27r)3 5(3) (k^ + + kg) [2 (ki, k^) V^{h)V^{h) + cyclic 

'.14) 

where V,i,{k) is the power-spectrum of the gravitational potential. 

At this point, in order to give the non-linearity parameter, an important remark is 
in order. Indeed, when dealing with second-order perturbations which are expressed 
in terms of first-order quantities, also the short-wavelength behaviour of the first-order 
perturbations must be taken into account, as it becomes evident when going to momentum 
space. The crucial point here is which is the final quantity one is interested in. We are 
interested in calculating the bispectrum of the gravitational potential on large scales as 
a measure of non-Gaussianity of the cosmological perturbations on those scales. The 
bispectrum of such quantities is twice the kernel which appears when expressing the 
second-order quantities in terms of first-order ones in Fourier space, an example of such a 
kernel being (ki, k2) in Eq. ()4.13|1 . This means that, when calculating the bispectrum, 
we can evaluate the kernel in the long- wavelength limit, irrespective of the integration over 
the whole range of momenta. This is the reason why the last term in Eq. ()4.12|1 gives a 
negligible contribution to the large-scale limit of the gravitational potential bispectrum. 
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Therefore, going to momentum space, from Eq. ()4.12|1 we directly read the corresponding 
non-hnearity parameter for scales entering the horizon during the matter-dominated stage 



r6 1 ki-k2 (ki-ka)^ ?>k'\ + ki , , 

fL - -2 + 4^ - 3^^^ + 2^fcr^ (4-15) 

where k = |ki + k2|. 

The non-Gaussianity provided by expression ()4.15|) will add to the known Newtonian 
and relativistic second-order contributions which are relevant on sub-horizon scales, such 
as the Rees-Sciama effect [18], whose detailed analysis has been given in Refs. [19]. 



5 Conclusions 

In this paper we have provided a framework to study the evolution of non-linearities 
present in the primordial cosmological perturbations seeded by inflation on super-horizon 
scales. The tiny non-Gaussianity generated during the inflationary epoch driven by a sin- 
gle scalar field gets enhanced in the post-inflationary stages giving rise to a non-negligible 
signature of non-linearity in the gravitational potentials. On the other hand, there are 
many physically motivated inflationary models which can easily accomodate for a pri- 
mordial value of /nl larger than unity. This is the case, for instance, of a large class 
of multi-field inflation models which leads to either non-Gaussian isocurvature perturba- 
tions [20] or cross-correlated non-Gaussian adiabatic and isocurvature modes [21]. Other 
interesting possibilities include the "curvaton" model, where the late time decay of a 
scalar field other than the inflaton induces curvature perturbations [22] , and the so-called 
"inhomogeneous reheating" mechanism where the curvature perturbations are generated 
by spatial variations of the inflaton decay rate [23]. Our findings indicate that a positive 
future detection of non-linearity in the CMB anisotropy pattern will not rule out single 
field models as responsible for seeding structure formation in our Universe. 
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